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Abstract 

The Riemann-Lanczos problem for 4-dimensional manifolds was dis- 
cussed by Bampi and Caviglia. Using exterior differential systems they 
showed that it was not an involutory differential system until a suitable 
prolongation was made. Here, we introduce the alternative Janet-Riquier 
theory and use it to consider the Riemann-Lanczos problem in 2 and 3 di- 
mensions. We find that in 2 dimensions, the Riemann-Lanczos problem is 
a differential system in involution. It depends on one arbitrary function 
of 2 independent variables when no differential gauge condition is imposed 
but on 2 arbitrary functions of one independent variable when the differ- 
ential gauge condition is imposed. For each of the two possible signatures 
we give the general solution in both instances to show that the occurrence 
of characteristic coordinates need not affect the result. In 3 dimensions, 
the Riemann-Lanczos problem is not in involution as a so-called "internal" 
identity occurs. This does not prevent the existence of singular solutions. 
A prolongation of this problem, where an integrability condition is added, 
leads to an involutory prolonged system and thereby generates non-singular 
solutions of the prolonged Riemann-Lanczos problem. We give a singular so- 
lution for the unprolonged Riemann-Lanczos problem for the 3-dimensional 
reduced Godel spacetime. 
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I. Introduction 



The problem of generating the spacetime Weyl conformal curvature tensor 
Cabcd from a tensor potential is called the Weyl-Lanczos problem and the 
analogous problem for the Riemann curvature tensor the Riemann-Lanczos 
problem. 

The work in this paper is based on the papers |||, ^ and on for 
the Riemann-Lanczos problem as an exterior differential system. But here, 
for simplicity, we mainly look at the Riemann-Lanczos problem in 2 and 
3 dimensions as an exterior differential system and then use Janet-Riquier 
theory to verify our results. As already explained in many papers such as 

|2|, 1^] we can express the Riemann tensor by means of a third-order tensor 
potential called the Lanczos tensor with components Late which are subject 
to the following symmetries 

L[abc] = (1) 

and 

Labc = L[ab]c ■ (2) 

If (|^) and (|2|) are imposed, we obtain 2 independent components in 2 di- 
mensions, 8 independent in 3 dimensions and 20 in 4 dimensions. The 
Riemann-Lanczos equations were first published by Udeschini Brinis Q as 

Rabcd = Labc-d ~ Labd;c + Lcda;b ~ Lcdb;a j (3) 

where ";" denotes covariant differentiation. However, as explained in 0], 
they did not form a system in involution but had to be prolongated once in 
order to be in involution. In what follows we include the differential gauge 
conditions 

Lab% = (4) 
but not the trace-free gauge conditions 

La% = . (5) 

If equations (^) were to be imposed, we would obtain the equation 

R = AL^\.n = -4L"\;fc = 

for the Ricci scalar R which leads to an inconsistent theory. 



II. Janet-Riquier Theory 

In the original Janet-Riquier theory of systems of partial differential equa- 
tions (PDEs), there was an algorithm created to explain how a given system 
of PDEs could be brought into passive form. Passivity is the absence of 
integrability conditions and orthonomicity is a form of ordering the partial 
derivatives of a system. The passive orthonomic system of PDEs was the 
predecessor of what is now called a formally integrable system or involutory 
system of PDEs. 

For a system of partial differential equations it is always an important 
question whether it possesses a general solution or whether further condi- 
tions have to be imposed so that the system is formally integrable or in 
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involution. A theory created by Riquier |22] and developed by Janet |11, 
and Thomas [13| now called Janet-Riquier theory helps to decide such ques- 
tions. Applications of the work of Riquier to general relativity have previ- 
ously been made in papers by Pereira |15, 
theory as in the form as created in 



17 p. A good account of this 



22, |ll|, 13] is given in Russian by S.P. 



Finikov |]. 

The original approach by Janet, Riquier and Thomas however did not 



lead to intrinsic results. To achieve this Spencer |28| and then Goldschmidt 
[10| introduced a new coordinate independent approach based on homolog- 
ical algebra. An account of this theory can be found in [19|. The earlier 
approach by Thomas [13| was thorough but based on elaborate systems of 
inequalities for monomials corresponding to a system of PDEs. 

There are many algorithms which have been implemented into algebraic 
computing. Some of them are concerned with the choice of particular rank- 



ings such as the REDUCE package DINV by Gerdt |2g]. Reid ^ developed 
a computer package based on MAPLE which brings a system into solved form 



which is a modification of an orthonomic system |20|. Seller [25, 26 1 uses the 
theory of the involutive symbol, which is a modernised version of the original 
Janet-Riquier theory, to determine whether systems are in involution. 

Next, we shall explain a few terms which are important in the original 
as well as in modernised versions of Janet-Riquier theory. For a system 
of partial differential equations (=system of PDEs) of order q, which we 
shall denote by TZg from now on, we use x^,---,x" for the independent 
variables, where we shall use brackets to indicate powers of any such as 
in (a;*)", and for the m dependent variables. Their derivatives on the 
jet bundle of order q, denoted by J7'''(M", M™), are denoted by u°j, where J 



[X 



"V" with 



is a multi- index. We associate a monomial x"^ = {x^ 
each partial derivative u°j, where J = {ji, ■ ■ ■ , jn) is a multi-index and 
ll'^ll = Ji + ■ ■ ■ + Jn is the order of the derivative. This means that to each 
set of partial derivatives of each dependent variable corresponds a unique 
set of monomials. 

Then, we order the partial derivatives of a system of PDEs in a sys- 
tematic way. Very often when x^, x'^, • ■ ■ , are our independent variables, 
we shall use an inverse lexicographic ordering based on x" ^ • • • ;^ x^ 1^ x^ 
on their partial derivatives so that y ■ ■ ■ y u^2 >~ u,i and so on. Gener- 
ally, we introduce a ranking amongst all partial derivatives 

Definition 1 Ranking of Derivatives 

A ranking of derivatives is a total ordering TZ-^ of all the partial derivatives 
u°j (for m,n fixed) satisfying the two conditions (where J and K are multi- 
indices) 

i) If\\J\\ > \\K\ 



then u°'j >- u 



a) If u' 



then 



K- 

K y {w 



),K for any multi-index K . 



A special subclass of these rankings are called orderly rankings by which 
we mean that these are rankings such that i) holds for different indices a 
referring to different unknowns u". A more detailed account on the problem 
of finding a suitable ranking especially for non-linear systems of equations is 
given by Rust [ 23 1 . Every ranking defined on the partial derivatives induces 
a ranking amongst the monomials which correspond to the given partial 
derivatives. 

Once a ranking for a system of equations TZq is determined, the system 
can be brought into a more organised form called orthonomic form. This 
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is achieved by determining the partial derivative highest in the ranking in 
each equation of TZq and calhng it the leading derivative of the equation. 
Once this equation is solved for its leading derivative which then becomes 
the only term on its LHS, we call it a principal derivative. All other partial 
derivatives of that order which are not in the set of principal derivatives are 
called parametric derivatives. Based on this we can define an 

Definition 2 Orthonomic System 

A system of partial differential equations TZq of order q is orthonomic with 
respect to a given ranking IZ-^ if 

i) all the PDEs are solved with respect to their leading derivatives; 
a) no two leading derivatives are the same; 

Hi ) and no parametric derivative in any equation of TZq can be a principal 
derivative in another equation of TZq or even a partial derivative of any order 
of a principal derivative. 

We base our calculations on definition (^) as our own standard. For or- 
thonomic systems the only term on each LHS is the principal derivative 
so that all derivatives of order q on the RHSs of the equations are para- 
metric derivatives. By differentiating each such equation we can also get 
the derivatives of the principal derivatives in such a form that we can de- 
fine multiplicative variables for each equation. Those independent vari- 
ables by means of which we can differentiate the principal derivative of an 
equation without re-introducing a derivative already produced by means of 
differentiating another principal derivative are called multiplicative vari- 
ables for the equation. In the same way, once we associate a monomial 
with each principal derivative, we can define what the multiplicative vari- 
ables are for that monomial. We say that a variable is multiplicative 
for the monomials which are of maximal degree in x". Further, x' with 
i < n is multiplicative for x"^ = {x^y^ ■ ■ ■ (x")-'" if amongst all monomi- 
als of the form (x^)'^i • • • (x*)'^*(x*"''-'^)-'*+i • • • (x")-'" the monomial x'^is such 
that ji = maxfc/cj. This definition given in p^] is based on Janet's original 



definition |1I]. 



We call a set of monomials complete if any multiple can be obtained 
using multiplicative variables only. We can also define completeness using 
partial derivatives. Given an orthonomic system, we denote the set of all 
leading derivatives and their derivatives with respect to multiplicative vari- 
ables only of the same dependent variable n" by . The closure then 
includes all derivatives of the derivatives of We call a system of partial 
differential equations of order q, TZq, a complete system when Cy = 
for all ti" , a = 1, • • • ,n . Once our system is complete, we can also ask 
whether it is passive. A system of PDEs having a complete set of monomi- 
als is called passive if any computation of a principal derivative is equivalent 
to any computation obtained using multiplicative variables only. Otherwise, 
the additional equations have to be added as integr ability conditions. If a 
complete system of PDEs is given in orthonomic form, we speak of a 

Definition 3 Passive Orthonomic System 

An orthonomic system of equations TZq is passive with respect to a ranking 
TZy, if the sets of all leading derivatives with respect to the chosen ranking 
TZ-^ are complete and no integrability conditions occur. 



Riquier |22f| formulated an important theorem for the local existence of an- 
alytic solutions based on this, namely, 
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Theorem 1 Riquier 

Given a passive system with respect to an orderly ranking TZy , then its formal 
power series solution converges. 

The main difficulty which arises in Janet-Riquier theory is that so many 
results are not intrinsic and depend on the choice of a ranking. In a more 
recent version this problem is re-examined by introducing special coordinates 
called (^-regular coordinates and by using the theory of the involutive 
symbol of a system IZq which we shall introduce next. But before that, 
we shall discuss a Janet example in detail in order to illustrate the theory 
exhibited so far. 

A classical example of Janet [12, 25, consists of two second-order 
equations with 3 independent variables and one dependent variable 

u so that n = 3, m = 1 and q = 2. The two equations are 

u,33 - x'^u,u = , 

U,22 = . (6) 

If we impose the coordinate ranking >- x^ >~ x^, we obtain the ranking 
"",33 >~ "",23 >- ^,22 >- "",13 ^ ^,12 >- ^,11 amongst the second-order partial 
derivatives. The leading derivatives then are 7/^33 in the first equation and 
u 22 in the second equation so that the system in orthonomic form is given 
by 

2 12 3 

U^33 = X U^ii XXX 

n 22 = X^X^ • . (7) 

The principal derivatives are ii^33 and u^22 which have the corresponding set 
of monomials {(x^)^, (x^)^} and the parametric derivatives are ii,i2i ^*,13) 
it^23- In the first equation all variables x^,x^,x^ are multiplicative variables 
but only x^ and x^ in the second equation whereas the variable x^ is not a 
multiplicative variable in the second equation (which is indicated by a dot). 
This is because, by means of differentiating the first equation twice with 
respect to x^, ti,2233 and higher derivatives can already be created^. For the 



Janet example (^, its completion is given in |19|,|£5|. When choosing the 



above ranking based on x^ ^ x^ ;^ x^, we obtain for jC^ and for jC^ 

'C>- = {'U,33;n,22;'U,133;'",233;'W,333;'W,122;'U,222} 
= {■U,33;n,22;'U,133;'W,233;'W,333;'W,122;'U,222;'U,223} • 

Therefore, ti,223 cannot be produced using multiplicative variables only - 
a fact reflected by the impossibility of producing the monomial (x^)^ • x^ 
using the set {(x^)^, (x^)^}. We have to add n^223 so that the new completed 
system is 

n,223 = x^ x^ • 

^,33 — X^U^ii =0 X"*^ x^ x^ 

?i^22 = X"*^ X^ • . 

But this system is not passive because, using multiplicative variables only, 
we can form n_2233 — 2;^n^ii22 — 2u^ii2 = which leads to the integrability 
condition n^ii2 = 0. We refer the reader to [19| for a detailed discussion of 
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the completion process and only state the final prolonged complete system 
as 



^^,11113 = 







• • 


W1123 = 







• • 


= 





x' 


• • 


H,223 = 





x^ 




U,112 = 





x' 


• • 


U,33 - x'^U^ii 


= 


x^ 


2 

X X 


U,22 = 





x^ 


X^ • 



where prolongation of a system of equations TZq simply means adding all 
derivatives of order g + 1 of the equations in TZq leading to the prolonged 
system of equations TZq+i, a process which can be carried out several times 
if necessary. 



A. Involutive Symbol and Formal Integrability 



In modern versions following Spencer ||2^ one wishes to replace the coordinate- 
dependent theory of complete passive orthonomic systems of equations TZq 
by means of quantities which will no longer depend on a choice of coor- 
dinates. A theoretical approach based on map diagrams from homological 



algebra |28, ^ also explained in [^] does fulfil these requirements. 

For practical purposes though we shall rely on the theory of the invo- 
lutive symbol and on 5-regular coordinates. The symbol of a system 
of PDEs involves only the highest-order partial derivatives of each equation 
in TZq which simplifies calculations especially of large systems significantly. 
We denote the quantities corresponding to each partial derivative ti'j by Vj 
and define the symbol Aiq as 

Definition 4 Symbol of TZq 

A system of partial differential equations TZq of order q locally described by p 
equations in solved form as ^'^ {x'' , , u°j) = forr = 1, has a solution 
space M.q for the unknowns Vf with a = 1, ...,m, \\J\\ = q: 

= 0, (8) 

where we formally differentiate with respect to the u°j. A4q is called the 
symbol of TZq. 

For simplicity the matrix rather than the map is usually regarded as the 
symbol of TZq. We associate with each symbol equation its multiplicative 
variables which are the same as the multiplicative variables its corresponding 
equation in TZq adopts. Once each symbol equation has its multiplicative 
variables allocated, we determine the class of an equation in A4q by 
counting the number of multiplicative variables it adopts - a number denoted 
by k such that < A; < n. We carry this out with all equations occurring in 
Mq, counting how many equations of each class there are. So we define 

Z?^'^^ := number of equations of class kinMq {oiTZq) . 
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The definition of the Cartan characters ^ Ug'^^is based on the Pg''^ and 
given by 

+ q — k — 1 



a 



1 



i''--=m-r^'-r']-Pi'\ (9) 



where m-Y^^^^ ("^g-i ^) the total number of partial derivatives of order q 

(k) 

that a system TZg will have. The Cartan character Oq represents the number 
of remaining independent partial derivatives of order q in the subsystem of 
class k after the removal of the number of the principal derivatives of class 

(k) 

k in that subsystem. In 5-regular coordinates aq equals the number of 
parametric derivatives of class k and order q. 

Because the notion "class of an equation" is ranking-dependent, we 
must ensure that we obtain intrinsic results. This can be achieved when 

(n) 

we use a ranking for which we obtain the maximal possible value for Pq 
and then for /Jg" + and so on until we obtain the maximal possible 
value for Ylk=i f^<t^ ■ This is equivalent to obtaining the minimal possible 
value for , then for a^" + and so on until the minimal possible 
value for Ylk=i ^^Q^ achieved. In an arbitrarily given coordinate system, 
a prolongation to higher order is sometimes necessary to obtain intrinsic 
results. This can be avoided by performing a coordinate transformation, 



where a linear coordinate transformation is sufficient |19], and by checking 



again for minimal and maximal values respectively of the a'q'^ and the (3^'^ 
in the above sense. Once this is fulfilled, we are using a system of (5-reguIar 
coordinates. 

It is important to determine whether a given system of equations TZq 
possesses identities or not. When no identities are present, the symbol is 
said to be involutive which can be equivalently expressed as 

Theorem 2 Involutive Symbol 

In a 5-regular coordinate system the following conditions are equivalent: 
i) The symbol A4q is involutive ; 

a) dim {Mq+i) = I]fc=i k ■ Oif^ ; 

Hi) for the rank r of A4q+i it is r{A4q^i) = J2k=i ^ ' (^'t'^ '■> 

iv) prolongation with respect to non- multiplicative variables does not 

lead to any new equations. 

But a system of equations TZq which has an involutive symbol Aiq can still 
admit integrability conditions. They can be revealed by means of projecting 
our prolonged system IZq+i, which is obtained by differentiating TZq with re- 
spect to all its n independent variables, back onto TZq. In general, we shall de- 
note first projections onto lower-order systems by 7^^+^ = 7rlj^_^J'~^^\TZq+r+i) ■ 

If TZ^g^ is not identical to TZq, then integrability conditions occur, but if they 
are identical we can characterise a formally integrable system TZq as ||2^, 



Definition 5 Formal Integrability 

A system of partial differential equations TZq is formally integrable means 
that TZ^q}_j. = TZq+r for all r > . 

A special situation occurs when Aiq is involutive and TZq is formally inte- 



grable [25|. Then, we obtain ideal systems of equations which are called 
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Definition 6 Involutive Systems of Equations 

A system of equations TZq is involutive if and only if its symbol Mq is invo- 
lutive and TZq is formally integrahle. 

Note that if we are dealing with a vector bundle, it is always valid dim(7^q"'^'') = 
dim(7^q+i) — dim(A^q+i). This means that when we know the dimension of 
the prolonged system of equations 7^g+i and the dimension of its prolonged 
symbol we automatically know the dimension of the system of equa- 

tions projected back onto order q called TZ^q \ This can be a great help in 
deciding whether a system has integrability conditions. 

The subsequent theorems, which can be found in are very useful 

for practical applications as only a limited number of steps need to be carried 
out to check for involutivity. First we look at 



Let Aiq be an involutive symbol of a system of equations TZ, 



Theorem 3 

then 

i) the symbol A^g+i is also involutive, and 



a) {TZq ^]^-^ — 'v-q+i 

Here i) states that, for any involutive symbol A4q, a prolongation is trivial 
and the prolonged symbol is again involutive and (ii) says that in this case 
projections and prolongations commute. A proposition following from this 
is 

Theorem 4 Proposition 

Given a system of differential equations TZq, where the number of equations 
equals the number of dependent variables m, then, TZq has either identities 
or integrability conditions which implies that (5^ > 0. 

Based on the previous definitions, we would have to perform infinitely many 
prolongations in order to decide whether a given system TZq is involutive or 
not. Fortunately this is not necessary and we can use the important 

Theorem 5 Criterion for Involution 

TZq is a system in involution means that its symbol A4q is involutive and 



7^' 



(1) 



TZ, 



(1) 



TZr, 



q — i^q- 

We shall now present the modern version of Janet-Riquier theory as 
given above using the Janet example (pi). Its symbol A42 is given by 



1^22 







1 2 3 

iJC tJO wjj 

1 2 
X X • 



where the integers on the RHS of each equation indicate the class of this 
equation. This leads to [3^^^ = 0, = (3f^ = 1 and so YL\^i k ■ (3^2'^ = 5 
which gives us the total number of multiplicative variables for M2 to 
be 5. In order for A42 to be involutive, we have to verify that r(Al3) is 
equal to 5. Firstly, we determine the new prolonged system of equations TZ^ 
which consists of the 8 equations 














■",33 - 

-",22 

"",133 

"",233 

■",333 

■",122 

■",222 

■",223 



X U iii 



X U^ii2 - U^ii 
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Next, we determine dim(7^3), where dim(7^g) is defined as the total number 
of dependent variables and all their partial derivatives up to order q minus 
the number of independent equations which constitute TZq. In this way, we 
can easily see that dim(7^3) = 20 — 8 = 12, where 8 is the total number of in- 
dependent equations in TZ-^. The number 20 stands for 20 = 1 + 3 + 6 + 10 for 
one unknown n, 3 first-order partial derivatives ti,2, ^^.Si 6 second-order 
partial derivatives it,ii, it,i2, ^*,13 5 "",22, ii,23, ^*,33 and 10 third-order partial 
derivatives n^m, n,ii2, n^ns, n,i22, m,i23, "",133, ""222, ^^223, "",233, ^i,333■ Next, we 
determine M-i which is given in orthonomic form as 



V333 


= x^Vm 




x^ 


x^ 


3 


V233 


= x'Vu2 




x"^ 


• 


2 


V223 


= 






• 


2 


V222 


= 


x^ 


x^ 


• 


2 


^133 


= x^Viii 


x^ 


• 


• 


1 


^122 


= 


x^ 


• 


• 


1, 



where we can see that there are 2 equations of class 1, 3 equations of class 2 
and 1 equation of class 3 resulting in = 2,(3^^ = 3,P^^ = 1. This leads 

to ^■/^s'^^ = 1-2+2-3+3-1 = 11 being the total number of multiplicative 

variables which Ai^ adopts. We also find that ^(Als) = 6 which is not equal 
the total number of multiplicative variables A42 adopts, which is 5, and 
therefore Ai^ is not involutive. We need to prolongate TZ^ further to TZ4 
and we refer the reader to |Q for further details subsequent calculations. 
We only state here that during this process the two integrability conditions 

-",112 = 

n 1111 = (10) 

occur. Finally, we obtain the result that the prolonged fifth-order system, 
which is augmented by (|lO|), 7^5^\ forms a system in involution and the 
Janet algorithm terminates. 



III. The Riemann-Lanczos Problem in 2 Dimen- 
sions 



We shall apply the theory of EDS and Janet-Riquier theory in order to 
examine the Riemann-Lanczos problems in 2 and in 3 dimensions. For a 
review of the theory of EDS see |, |, |2^, |0|, |, 0]. 

We show that the Riemann-Lanczos problem is in involution with respect 
to its two independent variables for both possible choices of signature and 
show that the occurrence of characteristic coordinates need not affect the 
result. We shall look at the problem first as an exterior differential system 
and then as a system of partial differential equations (systems of PDEs) 
applying Janet-Riquier theory. 



A. The Riemann-Lanczos Problem in 2 Dimensions as an 
EDS 

Firstly, we consider the Riemann-Lanczos problem in a 2-dimensional space- 
time. We know that in 2 spacetime dimensions, there is only one indepen- 
dent component of the Riemann tensor and therefore only one independent 
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Riemann-Lanczos equation. We now choose (x^, x^, L121, -L122) -Pi2ii; -Pi2i2) 
-Pi22i5 -P1222) as local coordinates on the jet bundle jr^(R^,]R^) so that it is 
8-dimensional. We leave the metric tensor components completely arbitrary 
so that the following results will hold for any 2-dimensional spacetime. The 
only independent Riemann-Lanczos equation is given in solved form as 



AR) 
J 1212 



-R1212 — 2P1212 + 2P1221 + 2Li2i(rj2 — r22) — 2Li22(rj;;^ + r^2) 
0. (Ill 



We define ai2i2e in the 2-dimensional case to be 



"12126 



R 



1212,e 



2Li2i(ri 



12,e 



22,e; 



2Li22(rh„ + r^ 



12,e; 



(12) 



where e = 1,2. The Pfaffian system derived from (O) is locally given by 



9^ = ai2i2idx^ + ai2i22C^x2 - 2dPi2i2 + 2dPi22i + 2{r\2 - r2 



22)dL 



121 



-2{T\, + Tl2)dLi22 , 



Pi2Mdx - P^2^2dx 



= dLi2i 

0^ = dLi22 - Pl22ldx^ - Pl222dx^ 

de"^ = dx^ A P1211 + dx^ A P1212 , 
de^ = dx^ A P1211 + dx^ A P1212 , 



(13) 



where 9^ is the exterior derivative of f^}^ and 6"^, 9^ are the two contact 
conditions and d9'^,d9^ their exterior derivatives. The zeroth character is 
So = 3, since we only have to count the number of 1-forms in (p!3|). Omitting 



all the terms involving dx^,dx'^ in the 1-forms in 13, the same number is 



obtained so that Sq = 3 leading to our first result sq = Sq = 3. In order to 
obtain si, we must form the first polar element H{{E^)x) of a 1-dimensional 
integral element {E^)x formed by the Vessiot vector field U 



U 
where 
because 



_d_ 



U'^ + PuieU' 



d 



dL 



+ Pl22eU' 



d 



121 



dL 



+ U- 



d 



121e 



122 



9Pi 



d 



122e 



121e 



f/1212 = C/1221 +(^1^/^ +<52f/^ 



9Pl22e 
(14) 

(15) 



61 
62 



«12121 + 0^12 ~ r22)-pL211 — {^U + ri2)Pl221 , 
"12122 + 0^12 — r22)Pl212 — {^11 + r^2)Pl222 • 



(16) 



This means that for the Vessiot vector field U the components U^,U'^, ^^1211, 
^^1221 , ^^1222 can be chosen arbitrarily and so we can for:mH{{E^)x). For the 



coefficient matrix of the second polar system H{(E'^ 


)x), we 


obtain 






equation 


dx^ 


dx'^ 


dLi21 


dLi22 


dPi2ii 


dPl212 


dPl22l 


dPl222 


"J 1212 

02 

2{U\d9^) 
2{U\d9'^) 
2{V\d9'^) 
2{V\d9^) 


"12121 
-P1211 
— P1221 

— f^l211 
— f^l221 

— V12II 

— ^1221 


"12122 

— P1212 

— P1222 
— f^l212 
— t^l222 

— V1212 

— ^1222 


'^{^12 ~ ^22) 
1 








-2(rli + r?2) 



1 



















-2 










2 




f/i 











y2 
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Here, U and V are the first and second Vessiot vector fields to be chosen. 
From the first 5 rows and their reduced counterparts of the above Matrix 
we determine si and s'^. We can easily see that the rank of this 5x8-Matrix 
is 5 and the rank for the reduced matrix, where we just omit the first 2 
columns for dx^ and dx"^, is 5 as well. Therefore, the Cartan character si, s'l 
are si = s'^ = 5 — 3 = 2. Apart from the same conditions that already 
hold for U, the second Vessiot vector field V must also satisfy the additional 
conditions 

V'Uabce - U'Vabce = . (17) 

For ^1212) the same conditions apply as for ?7i2i2 so that 

^1212 = ^1221 + SiV^ + hV^ (18) 

with 81,62 as defined above. The only free components left for V are V^,V'^, 
1/1221 • In order to obtain the rank of the second polar element, we need to 
consider all rows of the above matrix. First, we look at the reduced 6x7- 
matrix therefore omitting the first two columns under dx^ and dx'^. Let us 
denote the i^^ row by (i), i = 1, • • • , 7 and check whether the 7 rows are 
linearly dependent by looking for possible linear relations among them. We 
find 

(1) = Ai ■ (2) + A2 ■ (3) + Bi ■ (4) + B2 ■ (5) + Ci ■ (6) + C2 ■ (7) , (19) 
where 

A, = 2{T\,-rj,),A2 = -2ir\,-Ti,), 

Bi = 2y^, B2 = 2yi, 

L-i — > ^2 — ) 

and 6 = U'^V^ - U'^V^ . This leads to the result ^2=6 -3-2 = 1. Now, we 
must check whether or not the same linear combination occurs amongst the 
7 rows for the full polar system: 

(1) = Ai ■ (2) + A2 ■ (3) + Bi ■ (4) + B2 ■ (5) + Ci ■ (6) + C2 ■ (7) , (20) 

which must also hold for the first two columns. This leads to the two 
conditions 

ai2i2i = —Ai ■ P1211 — A2 ■ /'mi — Bi ■ i7i2ii — B2 ■ ?7i22i — Ci • V1211 
—C2 ■ V1221 ) 

ai2i22 = —Ai ■ P1212 — A2 ■ -P1222 — Bi ■ C/1212 — B2 ■ C/1222 — Ci ■ V1212 

—C2 ■ V1222 ) (21) 

where we have to insert the above known expressions for ?7i2i2, ^1212, ^1211 
and V1222 and it is not permitted to restrict any of the U^^V^. By inserting 
all the above determined values for the Ai,Bi,Ci, we see that in the linear 
relations (^T|) the free components ?7i2ii, C/1221, ^^1222, V1221 cancel out and 
both relations in (|2l| ) hold identically without imposing any further restric- 
tions. Therefore, the same linear combination ( [20| ) holds for the full system. 
The second character and its reduced counterpart therefore coincide so that 
S2 = 82 = 1 and the involutive genus g is g = X]i=i = 2. 

We can also compute the characters using the tableau and by determin- 
ing the torsion we then find out whether the system possesses integrability 
conditions. Given a Pfaffian system V consisting of s 1-forms 0° with in- 
dependence condition Q = uj^ A ■ ■ ■ A u>^ , we denote by vr^ all the extra 
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forms such that (^", w*, vr'^), 1 < a < s, form a coframe on our formahy 
A^-dimensional manifold Ai. We can then write 

= A'^i aJ + A w^' + ^C^^vr^ A t:^ (mod T{V)) . (22) 



In equations (p2D, the A"- form the tableau matrix and the -Bj" are caUed 
the torsion terms. Note that if the coefficients C"^ = 0, then the sys- 
tem is said to be quasi-linear. In order to form a complete coframe for 
our example (|l^), we see that we have to add three 1-forms vr^, where A 
now is a collective index with A € {1211,1221,1222}, to the five 1-forms 
9^ , 6'^ , 9^ , uj^ , u}"^ and so we choose 

7ri2ii := dPi2ii, 7ri22i := dPi22\, 7ri222 := d-Pi222 • 
The dPi2i2 can be expressed through 0^ as 

1 1 

dPi2l2 = ^1221 + (2ai2121 + (ris - G^2)(A211 - + Tl2)Pl22lW + 

ai2122 + {r\2 - G^2)(^1212 " {^li + rl^)Pl222)iO^ (mod {^j) (23) 

Having accomplished this, we can then write the exterior derivatives of the 
contact conditions as 

d9^ = 

d9'^ = -7ri2ii Auj^ - 7ri22i Auj^ + (-012121 + (r}2 - r22)-Pi2ii 
-(rn + r?2)Pi22i)cu^ Au;2 

d9^ = -7ri22i Alo^ - 7ri222 A . (24) 
Therefore, we obtain for the tableau matrices 



Aa 
^Al 



where a = 1, 2, 3 and A is one of the collective indices A G {1211, 1221, 1222}. 
This leads to s'^^ = 2, §2 = 1, and, the only non- vanishing torsion term is 
given by 

Bu = - [^«i2i2i - (r}2 - r^2) A211 + (rJi + r22)Pi22i] • (25) 

If we wish to absorb the torsion coefficients, we must find a transformation 
$ with vr'** — > vr'^ + p'\uj^ and quantities pf such that 

In our case this leads to the system 

= ^A2P^-^AlP^ 

= Al^p\-Al,p^^. (26) 

One solution Ansatz to fulfil (^) is to choose p\ := Bf2 while all other p^ 
vanish: thus the torsion is absorbed and the system is therefore in involution 





00^ 













00 


A12 = 





-1 




\ 


-10) 
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Adding the exterior derivative of the differential gauge condition 
Li2'^.,, = 0, which is locahy given by 

dLi2% = g^^dPun + g^HdPi22i + dPuu) + g'^^dPi222 - SadL^i - 5idLi22 
+5^dx^ + d^dx^ , (27) 

where S3,---,6q are defined in Appendix A, does not change the results 
qualitatively. We can write down the polar systems again and from this 
obtain the polar matrices of which the matrix of H{E'^)x is given in Appendix 
A. 

Obviously sq = Sq = 4 and we find that no linear combinations are 
possible amongst the first 6 rows or their reduced counterparts so that 
si = s[ = 2. But for S2 there exist multipliers Ai,Bj such that 

(7) = ^1 • (1) + A2 ■ (2) + ^3 • (3) + ^4 • (4) + A5 ■ (5) + • (6) , 

(8) = Bi ■ (1) + B2 ■ (2) + B3 ■ (3) + ^4 • (4) + B5 ■ (5) + Be ■ (6) (28) 

now with multipliers Ai, ■ ■ ■ , Aq and Bi, ■ ■ ■ , Bq as given in Appendix A. In 
order for these linear combinations to hold for the full rows (7) and (8) as 
well, the following 4 equations have to hold: 

14211 = -A\aY2\2\ - ^2(55 + ^3-Pl211 + ^4-Pl221 + A^\J\2\\ + A^ 

f/1221, (29) 

y\22\ = —-81012121 - -82^5 + -B3P1211 + B4P122I + -85^^1211 + Bq 

t/1221, (30) 

^1212 = — ^1^12122 — ^2^6 + ^3-^1212 + ^4-Pl222 + A5U1212 + Aq 

?7l222, (31) 

^1222 = —-81012122 — -82^6 + -B3P1212 + B4P1222 + -85^1212 + Bq 

«7l222, (32) 

where U is the first Vessiot vector field and V the second one. After a 
lengthy calculation of which details are given in Q we find that with ([l8| ) 
the above expressions (^) to ( |3^ ) hold. Therefore, the full rows (7) and (8) 
are again linearly dependent which leads to S2 = ^2 = and g = 2. Note 
that all the above results can be verified using a REDUCE code given in M. 



B. The 2-dimensional Riemann-Lanczos Problem as a Sys- 
tem of PDEs 

The Riemann-Lanczos problem in 2 dimensions is an example of a system 
of partial differential equations for which applying Janet-Riquier theory is 
much more economical. We shall also see in Section C below that it is 
of great importance to compute the characters using (5-regular coordinates. 
There, if characteristic coordinates are used, they will lead to non-intrinsic 
results so that a coordinate transformation will then have to be carried out. 
We have a system of equations TZi with n = m = 2,q = 1 and the symbol 
M.I consists only of a single equation 

- 2^1212 + 2^1221 = , (33) 

which is of class 2 whatever the choice of the orderly ranking. Note that 
here and in subsequent sections, where Janet-Riquier theory is applied, the 
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symbol variables Vabcd^yabcde have no connection to Vessiot vector fields V. 
This yields 

a« = 2, aP) = l. 

Trivially, all variables are multiplicative variables and the symbol of one 
equation is always involutive. But, we must examine whether this system 
is formally integrable. We only need to compute 7^2 and its projection TZ^^^ 
and check whether 7^2 = T^i^'') where dim(7^i) = 5 and dim(A4i) = 3. Our 
prolonged system of equations 7^2 consists of /i^i2 ' /1212 1 ' /1212 2 ' namely, 

= -R1212 — 2P1212 + 2P1221 + 2Li2i(rj2 — r22) — 2Li22(rj;i + rfg) 

= ^1212,1 — 2S'l2121 + 25i2211 + 2Pl21l(rj2 — r22) 

-2Pl22l(rii + r^g) + 2-^^121 (ri2 - ^22),! - 2Li22{T\^ + T\2),1 
= -Rl212,2 — 2S'i2122 + 2512212 + 2^1212(1^12 ~ ^22) 

— 2Pi222(rii + 1^12) + 2^121 (ri2 ~ 1^22), 2 — 2Li22(r|]^ + rf2),2 • (34) 

No integr ability conditions can be created because 

dim(7^^^^) = dim(7^2) - dim(A^2) = 9-4 = 5 = dim(7^l) 

and the system is linear so that formal integrability has been shown and the 
system TZi is involutive. 

If we incorporate the differential gauge condition Lab^.^ = 0, then M.i is 
formed by 

- 2V1212 + 2^1221 = 
5''^1211 + 5'Vi221 + 5^^1212 + 5'Vi222 = , (35) 



which has two equations of class 2 again whatever the choice of the orderly 
ranking be and we obtain = 2, ) = . In order to check whether M.i 
is involutive, we prolong to the corresponding M.2 which is formed by 



— 2^12112 + 2T/i2211 — 



-2142122 + 2Vi2212 
>112 

AUr_ _ , „12t7__ _ , „12t/_ , „22i 



5"'^^12111 + 5^^^12211 + 5^^Vi2112 + 5^Vi2212 = 



V12112 + 5 V^12212 + 5 V^12122 + 5^12222 = , (36) 

We get r(A^2) = 4 = A: • /jj'^'* and therefore Adi is involutive. We must 
check that the new system TZi is formally integrable again and so add the 
following 3 equations Li2*s, 9il/i2*s, 52^12*^5 to (^): 

= 5"'^^-Pl211 + 9^'^{Pl221 + -P1212) + 5^^-^1222 — Li2l5z — -^122(54 

= 5'"'^^'S'l2111 + 5"'^^('S'l2211 + 'S'12112) + g^'^Si22l2 + 5^1-^1211 + ^^^l (-^1221 + -P1212) 

+5^1-^1222 — -Pl211'53 — -Pl221<^4 — -^121(53,1 — ^122(^4,1 
= 5"'^^'S'l2112 + 5"'^^('S'l2212 + 5'i2122) + 9^"^ S 12222 + 5'^2-f'l211 + 9^'^{Pl221 + A212) 

+5^2-^1222 — -Pl212'53 — -Pl22254 — -^121(^3,2 — ^^122^4,2 , (37) 



where 5-\ , 6a defined in Appendix A. The system (p4) together with 

( ^7|) now forms the prolonged system of PDEs 7^2 including the differential 
gauge condition and its partial derivatives. Again, we cannot solve any of 
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the 4 second-order PDEs in such a way that no second-order derivatives 
remain. In order to show this, we can use the symbol matrix of as welh 



^12111 


^12112 


^12122 


^12211 


^12212 


^12222 





-2 





2 














-2 





2 





5" 


9'' 





9'' 


9'' 








9'' 


9'' 





9'' 


9'' 



which has maximal rank r(A42) = 4 and no integrability conditions can 
arise because of the linearity of the system. This means that TZ^^^ = 7^l, 
thus making TZi formally integrable and involutive. We conclude: 

Theorem 6 Proposition 

In any two-dimensional spacetime, the Riemann-Lanczos equations are in 
involution with respect to the spacetime coordinates and the differential gauge 
condition does not affect the result. The corresponding Cartan characters 
(sO) -51,52) for the involutory systems are (3,2,1) when no differential gauge 
condition is imposed but (4,2,0) when the differential gauge condition is 
imposed. 

C. General Solution to the Riemann-Lanczos Problem in 2 
Dimensions and Characteristic Coordinates 

It is instructive to repeat some of the above calculations using different local 
coordinates. First, we look at the general 2-dimensional spacetime with 
Lorentzian signature and line element in characteristic coordinates 
written as 

ds^ = -e'^Pdx^dx'^ , 

where p is an arbitrary function of x^,x'^. The only two equations we have 
to solve are /1212 — 0' -^12*5 = : 

= -e^Pp^^i^2 + 4Li2ip,^2 - 4Li22P,ai - 2P1212 + 2P1221 , 

= 2e-2^(2Li2i/9,^2 + 2Li22/0,^i - P1212 - A221) • (38) 

We find that the general solution is of the form 

Li2i = e2^/i(:ci)-ie2>,.i, 

L122 = e^'f2{x^) + -/'P,.^ , (39) 

where /i(x^), /2(a;^) are 2 arbitrary functions depending on one local coor- 
dinate each. This is in agreement with the result for the Cartan characters 
claiming that the general solution depends on two arbitrary functions of one 
variable each. But can we obtain the intrinsic values for the characters using 
the above coordinate frame? The symbol derived from ( |38|) is 

= -2^1212 + 2F122I 

= -1^1212 - ■^1221 , (40) 

from which we obtain (3^^ = (3^^ = 1 based on the ranking stemming from 
x"^ ^ x^. This leads to a[^^ = a^^^ = 1, which cannot be the intrinsic result 
as we already solved the 2-dimensional problem in the previous section. We 
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see that in ( |4(]| ) only one variable of class 2, namely V1212) and only one 
variable of class 1, namely 1^1221, occurs whereas neither ¥1211 nor V1222 
occur at all to make both equations be of class 2. As this can happen in 
another coordinate system as we shall see for the 3-dimensional case below 
as well, the above coordinate frame is not 5-regular and we need to perform 
a coordinate transformation of the form 

dx^ = a2idx^ + a22dx^ (41) 

in order to obtain the correct values for , af'^ . After such a transforma- 
tion, the new symbol in orthonomic form is 

^1222 = ^^1221 

022 

^1212 = -—^1211 . (42) 
ai2 

Now, both equations are of class 2 so that = 0, f3^^ = 2 which pro- 
duces = which is the known intrinsic result for the Cartan 
characters. 

For spaces with Euclidean signature we can write their line element as 

ds^ = e'^P{{dx^f + {dx^f) , 

where p is again an arbitrary function of x^,x'^ leading to the system of 
equations 

= -e^''(/3 + p^a:^x2) + 4Li21/9 ,^.2 - 'iLi22P^x^ - 2Pi212 + 2Pl221 , 

= e-^P{-2LuiP,xi - 2L122/0 .2 + P1211 + P1222) • (43) 



The solution to (43) resembles the one for (|3q) and is 



L,2i = e^'^{h{x' + x^)-^p^,2) 

L,22 = e'P{f2{x'-x^) + ^p^,^). (44) 

In this case we find the intrinsic values for the characters to be a[^^ = 2, 
af'^ = using ( [43| ) directly due to the fact that here we did not use charac- 
teristic local coordinates. 



IV. The Riemann-Lanczos Problem in 3 Dimen- 
sions: Involution or Prolongation? 

In this section we shall discuss the Riemann-Lanczos problem in 3 dimen- 
sions as an EDS and as a Janet-Riquier system of partial differential equa- 
tions (PDEs). First, we shall use EDS theory to find out more about the 
problem. 

A. The 3-dimensional Riemann-Lanczos Problem as an EDS 

In 3 spacetime dimensions with local coordinates x^, x^, x^, we obtain 8 in- 
dependent components of the Lanczos tensor, namely, L121, L122, Li^i, L133, 
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-^232) -^233 5 -^123 7-^132 when imposing the cydic conditions (||). Each of them 
has three first-order partial derivatives so that we use the jet bundle 
with formal dimension = 35 to express the EDS. We obtain the 6 indepen- 
dent Riemann-Lanczos equations whose exterior derivatives are df [212 ; 4f 1313 > 
'^/2323 ' ^fuis ' ^/i223 ' ^/i323 • These latter equations can locally be written as 

= -2dPi2i2 + 2dPi22i + dLi2ii2T\2 + 2TI2) - dLi22i2T\^ + 2r?2) + 2ri2 

■dLisi + 2r?idL232 + 2rf2C^^123 " 4r?2(iLi32 + (i?1212,e + 

= -2dPi3i3 + 2dPi33i + dii3i(2rl3 + 2ri3) - dL^si^Tl, + r?3) + 2ri3 

■dLi2i - 2r?idL233 - 4Tl^dLi23 + 2r?3dLi32 + (i?1313,e + ai313e)rfx^ 
= -2dP2323 + 2dP2332 + 2dL232(ri3 + T^g) - 2(iL233(ri2 + ^Is) " '^^h 

■dLi22 - 2Tl2dLr33 + 2r23dLi23 + 2Tlr^dLi32 + (^2323,e + a2323e)c^a;^ 

= -dPi213 + dPi231 - dPi3i2 + dPl321 + dLi2l(r}3 + 2ri3) + dLisi 

■(^12 + 2r23) — rf2d-^^133 — r^3(iLi22 — r^]^(iL232 + TlidL233 — dLi23 
•(rJi + 2r22 - r^g) - dLi32(2rf3 + - rj^) + (iil213,e + ai213e)rfx^ 

= -dPi223 + 2fiPl232 - C^A322 + ^^^2321 + f^il22(2rl3 + T^g) - dL232 
•(^12 + 2TI^) - Tl^dLi2l + rf 2^^233 + r22^131 " r22C?-^^133 " dLi23 

•(rL + 2ri2 + ri3) + ciLi32(2r33 + 

^22 ~ ^^2) + (-Rl223,e + Oli223e)dx^ 
= -2(iPi323 + dPis32 + dPi233 + dP2331 " 'i^l33(2ri2 + ri3) - dL233 

'{2^12 + ^13) + Tl^dL232 + r|3(iLi22 + r23Li3i — Tl'^dLi2l — dLi23{2Tl^ 
-ri3 - ris) + dLi32{r\s + 2TIs + r23) + (i?i323,e + ai323e)(ix^ , (45) 

where aabcde is defined in Appendix B. The polar matrix consists of 35 
columns for the {dx^ ,dLabc,dPabcd) and we split the first 6 rows stemming 
from (p5|) into 3 blocks which we call (Mqi Mq2 Mq-^) so that the full 
polar matrix of H{E^)x is given by 

/Moi M02M03 \ 
M(P) 1 
M{U) N{U) 
M{V) N{V) 
\M{Z) N{Z) J 

where we denote the truncated matrix obtained by leaving out the rows 
(Mqi Mq2 Mq^) by Mi. All the matrix expressions can be found in Appendix 
B and a more detailed explanation is given in The first 8 rows from Mi 
together with (Mqi -M02 M03) lead to sq = s'q = 14. Then, for si and s'l 
there is neither a linear combination of the full polar system nor one of the 
reduced polar system possible. We have si = s'l = 8 because the rank of the 
relevant 22x35-matrix consisting of (Mqi Mq2 M03) and the first 16 rows 
of Ml is 22 = So + si = Sq + s'l and so sq = s'q = 8. But now, for the second 
reduced polar system, we find that there is exactly one linear combination 
possible which we can write down as 

= ^1 • (rowi) + A2 • {row2) + ■ ■ ■ + Aq ■ {rowQ) 

Bi ■ (rowr) + B2 ■ (rows) H \- Bs ■ (rowu) 

+Ci ■ (rowi^) + C2 ■ {rowio) H 1- C's • {row22) , (46) 

where the linear multipliers are given in Appendix B. If we wish that this 
linear combination for the full polar system be valid, we have to claim that 
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the 3 equations (for e = 1, 2, 3) hold 



= ai212e^l + ai313e^2 + "23236^3 + "12136^4 + Q;i223e^5 + "13236^017) 

— f^l21e-Bi — Ui22eB2 — C^131e-B3 — C^133e-B4 — ^72326-^5 — f/233e-B6 

— Ul23eB7 — UiS2eB8 — Vi21eC'i — ^1226^2 — ViaigCs — 14336^^4 

— ^2326^5 — V233eC6 — ^1236^7 — Vi32eC8 , 

where Oabcde is as defined in Appendix B, without imposing any further 
restrictions on the first two Vessiot vector fields U, V. This can easily be 
checked by looking at how many of the Uabcd and how many of the Vabcd can 
be chosen arbitrarily. For the 24 Uabcd there are 6 restrictions, so 18 can 
be chosen arbitrarily. For the Vabcd we have 14 restrictions, so 10 remain 
arbitrary. Based on this argument, we find that some of the Uabcd, Vabcd never 
cancel out and the linear combination (^) does not hold for the full system 
so that S2 = S2 + I such that S2 = 7 and S2 = 8 |^] . This system is clearly not 
in involution and its reduced characters (sq,s^, 4,4) are (14,8,7,3). We 
can test this result with the REDUCE code given in |^] and immediately 
obtain (8,7,3) for the reduced characters {s'^,S2,s'^) and that the system is 
not in involution. In order to carry out this computation it is sufficient to 
use an arbitrary diagonalized line element which describes 3 -dimensional 
spacetimes completely as shown in 

We can also use the tableau matrix to calculate this result as shown 
next. First, we must complete our set of 14 one- forms together with the 3 
cj'^ to a complete coframe consisting of = 35 elements by introducing the 
18 vr^ such that {dfj^^^^,Kabc,^'^,'^^) forms a cobasis on J\4. We introduce 
the 18 vr"^ using collective indices A with A = 1, • ■ ■ , 18 with the following 
correspondence between vr^ ^ dPabcd'- 

TT^ ^ dPl21l , TT^ w dPi221 , VT^ ^ ^Pisn , 

vr^ ^ d-Pi33i , ^ dP2mi , Tr*^ ^ d-P233i , 

vr'^ ^ dPi23i , vr^ ^ dPi32i , vr^ ^ d-Pi222 , 
vr^o ^ dPi3i2 , vr^i ^ ^^1332 , vr^^ ^ ^^2322 , 

^1322 , 
^1233 • 



^ dP2332 , vri4 ^ dPi232 , vr^^ ^ dPr. 

7^16 ^ c?Pi333 , Vrl7 ^ dP2333 , VT^S ^ dPr. 



When using the correspondence ^" ^ Kabc based on equations (18) in 
the 8 dKabc can be recast as 



de"" 








- vr^ A - 


vr'^ A + vr^o A - vr^ 


AuJ^ 








AtJ^ H 




uj"- Au)'^ 




de^ 








-tt'Au^- 


27ri^Aa;3-vr5A^3 - 

2 


Bf3^ 


dB^ 






AuJ^ 


— TT 


Au? - 


- vr"^ A + ^Sfgu;* A tJ^ 




de^ 






Aoj^ 


— vr 


AlJ - 


- tt'^^ a uj^ 




dO^ 






Auj^ 


— vr 


AlJ - 


- vr^^ Auj^ + ^B%uj' A oj^ 




dO^ 






AuJ^ 


— TT 


Au? - 


-tt^'Au;' 




de'^ 










Au? - 


-n'^Au;' 


-V 
2 


de^ 










AlJ - 


- -TT^l A - -TT^ AW^ - 

2 2 






+ 2^i3^ 


AuJ^ . 









(48) 



18 



We can easily calculate the reduced characters when using the tableau ma- 
trices derived from ( [48| ) and obtain s[ = 8, S2 = 7 and S3 = 3. We also find 



that all remaining torsion coefficients B'j^j in ( [48D can be absorbed so that 
no integrability conditions can occur for the system (^5|) itself. However, 
because we have the rank deficiency in the reduced polar matrix leading to 
S2 = ^2 + 1, we must prolongate the system to second order. We introduce 
new jet variables Sabcde, of which the projection back onto the spacetime 
manifold corresponds to Sabcde ■= Pabcd,e = Labc,de, and denote the 24 con- 
tact conditions arising from them by Kabcd- The prolonged Pfaffian system 
then is given by 

^ ~ J abed 
^ ~ Jabcd,e 

= rf/S 

= df^Sd,e 
Kabc = dLabc ~ Pabcddx'^ 
Kabcd = dPabcd — Sabcdedx^ (49) 

together with the exterior derivatives dKabc and dKabcd- Calculations using 
an adapted REDUCE computer code suggest that the system (^9|) is not in 
involution and integrability conditions occur. 

Note that instead of the prolongation to ( ^9|) we can carry out the same 
prolongation which was carried out for the Riemann-Lanczos problem in 
4 dimensions in This may lead an involutory system eventually but 
we leave this as future work. However, the process of prolongation often 
becomes simpler when Janet-Riquier theory is used what we are going to 
next. 

B. The 3-dimensional Riemann-Lanczos Problem as a Sys- 
tem of PDEs 

The symbol ^Al for the unprolonged Riemann-Lanczos problem in 3 dimen- 
sions consists of the 6 linear equations derived from f[2i2 > /1313 > /2323 ' /1213 ' 
/1223 ' /i323 given in orthonomic form as 



V1212 


= ^1221 




1 2 
X X • 


^1313 


= ^1331 




1 2 3 

tJO 


V2323 


= ^2332 




1 2 3 

CO CO 'JO 


V1213 


= ^1231 — 


^1312 + ^1321 


1 2 3 

CO CO CO 


^1223 


= 2^1232 - 


- ^1322 + V232I 


1 2 3 

CO CO CO 


^1323 


= ^(^1332 


+ ^1233 + V2331) 


12 3 

CO CO CO 



We imposed the ranking x^ >- x"^ >- which then induced the ranking 
PabcS, >- Pabc2 ^ Pabci and amongst each set Pabce the ranking P233e ^ 

^^2326 >- P\32e ^ -pL23e ^ -Pl33e ^ Pl31e >- Pl22e ^ Pl21e- From this We nOW 

obtain (3[^' = 0,/3]^' = l,/3f^ = 5 and a^^ = 8,aY' = 7,af^ = 3. But, we 
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must check whether we actually used (5-regular coordinates and we perform 
the following coordinate transformation 



dx^ = andx^ + a^dx'^ + ai^dx^ 
dx^ = a2idx^ + a22dx'^ + a23dx^ 

dx^ = a-iidx^ + a32(ix^ + a^^dx^ , (50) 
which leaves us with the symbol in orthonomic form in the new coordinates 



^1223 — {Bi + 023^1213) 

ai3 

^1333 = {B2 + 033^1313) 

ai3 

^2333 = (-B3 H — ^(-85 H — —Bi H —B4^ 

023 «13 «13 «13 

+2V1213— ^^-^ — 3023^1233 H ^^1313) 

«13 0-13 

V1323 = {B4 + 033^1213 — «13^1233 + 023^1313) 

013 

V^2332 = /(Kfecl,^afec2). (51) 



In (^) we dropped the and the Bi and f{Vabci, V^6c2) are given in Appendix 
C. We see from (|5l| ) and Appendix C that no variable of the form Vabc3 
remains in the last equation for any choice of the aij in ( [50| ) so that we 
obtain 5 equations of class 3 and one of class 2. This means that our original 
coordinates were 5-regular and we can use our first coordinate system in 
order to obtain intrinsic results. 

In order to see whether the symbol is involutive, we differentiate each 
equation with respect to from which we can produce A42, where 

its sparse coefficient matrix is given and explained in detail in |p. We can 
determine the rank of this sparse matrix easily by hand (or using Maple) 
and obtain r{A42) = 18 as shown in We now compute the total number 
of multiplicative variables 



fc=i 



= 17 / r{M2) = 18 



which means that even the necessary condition for the symbol to be involu- 
tive is not fulfilled and the system cannot be in involution. 

Because A4i is not involutive, we prolong TZi to 7^2 which consists of 24 
equations 18 of which are the partial derivatives of the 6 Z^^^,^- We write the 
prolonged symbol A42 down which is based on the following ranking. We 

order the Sabcde such that Sabc33 >- Sabc23 >- Sabc22 >- Sabcl3 ^ SabcU >- SabcU 

and then amongst each set 52331^ ^ 5'232ii ^ Si32ij >- Si23ij >- Sr33ij >- 
Si3iij >~ Si22ij >- Si2iij- We find that the symbol in orthonomic form is 
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given by (with multiplicative variables indicated in each equation) 



T 


^12112 


= ^12211 






• 


• 




^12122 


= ^12212 








• 


"3 


^12123 


= 2V12312 


— ^13212 + V232II 


X^ 




• 



4 


V13113 


= ^13311 


X^ • • 




V13123 


= ^13312 


1 2 
X X • 


6" 


^13133 


= ^13313 


12 3 
X X X 


7 


V23213 


= V23312 


X^ • • 


8 


^23223 


= ^23322 


1 2 
X X • 


9" 


V23233 


= ^23323 


X^ X^ 





^12113 


IT 


^13122 


T2" 


^12133 



^12311 - 

2^13212 

1^12313 



^13112 + ^13211 
- ^12312 — V232II 



^13312 + ^^23311 



1 2 
X X 



1 9 3 

tJC wJO wJO 



13 


^12213 


14 


^12223 


15 


^12233 




^13213 


il 


^13223 


28_ 


^13233 



2V12312 — ^13212 + V232II 
2V12322 — ^13222 + V23212 
1^12323 + 5^23312 — ^^^13322 



17 ^13223 — ^(^13322 + ^^12323 + V23312) 



X • • 

1 2 
X X • 

12 3 
X X X 



X • • 

1 2 
X X • 

1 2 3 
X X X . 



This system produces 6 equations of class 1, 7 of class 2 and 5 of class 3 
leading to /Jg^^ = 6, /Jg^^ = 7, fi'^^ = 5. The total number of multiplicative 
variables equals 35 and if M2 were involutive, we would have to obtain 
^(-^3) = 35. Differentiating IZ2 with respect to all 3 spacetime coordinates, 
we obtain 54 equations relevant for the symbol But all those equations 
obtained from differentiation with respect to non-multiplicative variables 
are superfluous as each of them can be obtained by differentiating with 
respect to multiplicative variables except for c?3/{2i2- Therefore, r{Ais) < 
36, when we use formal differentiation of the above symbol equations directly 
because d3{eqn.3) corresponding to 03/^^12 is a linear combination of some 
other equations in ^^3. We obtain r(A^3) = 35, which means that A^2 is 
involutive. 

We obtain that the symbol equation d3{eqn.3) can also be created as a 
linear combination of 



d3{eqn.3) = 2d2{eqn.l6) + 2ds{eqn.ll) — 2ds{eqn.l3) — d2{eqn.5) — di{eqn.7) 

(52) 

using formal differentiation so that r{M.3) = Yli^k=i ^ ' f^i^^ ~ 35- However, 
it turns out that when (^) is rewritten in terms of the full equations from 
7^3 as 

j_ r{R) , fiR) I f{R) _Of(^) _9f(^) _l_9f(^) (t^'i'\ 

J 1212,33 ^ 1313,22 "T ./ 2323,11 ^./1323,12 1213,23 1223,13 ' W-JJ 

it is not a trivial identity any longer and cannot be obtained by means 
of any linear combination of the /^^|^ and their derivatives g so that 
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r{TZ^) = 6 + 18 + 36 = 60. Because our system is linear, it is instructive to 
compute 

dim(7ei) = 18 , dim(A^i) = 18 , 
dim(7^2) = 56 , dim(Al2) = 30 , 
dim(7^3) = 100 , dimj^Ig) = 45 , 

which leads us to the important consequence that 

dim(7^^^^) = dim(7e3) - dim(A43) = 100 - 45 = 55 = dim(7^2) - 1 . 



But this means that one integrability condition of the form (|53D oc- 
curs. Our new system of equations has to be given by 

~ J abed 

- f^^^ 
^ ~ Jabed,e 

= 1 . (54) 

We find that its symbol is simply given by the previous together with 

^12323 = p3 _-p3 [V23312(r?3-rf2)+2r|3¥L2312-^13322(2rJi+rf2 + rf3)] 
^12 ^13 

(55) 

which adopts 3 multiplicative variables. Therefore: 

f^Ui) = 6, /32^(\) = 7, = 6 so that a^^J-^^ = 18, afj^^^ = 9, a^^(\^ = 2 . 

We see that r{Ml^^) = r^M-s) + 3 = 38 because the 3 derivatives resulting 
from ( |55|) are linearly independent and because 

r(7^^^)) = r(7^l) + r(7e5^^) + 35 + 3 = 63 

which leads to 

dim(7^f ^) = dim(7e^^^) - dim(7W^^^) = 97 - 42 = 55 = dim(7^^^^) . 

Because r(X^^^) = ELi ^'/^S^) = 6 + 14 + 18 = 38, M^2^ is involutive and 
because dim(7^2^^) = dim(7^2^^) = 55 , there are no integrability conditions 
that can occur and so the system TZ'^^ is in involution. 

C. A Covariant Formulation of the Integrability Condition 

for n = 3 

Even though the above prolongation of the Riemann-Lanczos problem math- 
ematically leads to a prolonged system in involution, this prolonged system 
( |5^ does not respect general covariance. This is not satisfactory from a 
general relativity point of view and we are going to suggest a covariant ver- 
sion of the above prolongation. Instead of adding the partial derivatives to 
our equations, we look at the following new system of equations 

^ ~ J abed 

= /S;e ■ (56) 

The symbol of this modified system obviously coincides with that of the 
system formed by fj^f^^, f^}^ g and again, we obtain for the rank of ^^3 that 



22 



riMs) = Ylk=i k- = 35 so that M2 is involutive. We now have to see 



whether the hnear combination (52) also holds for the full covariant system 
here. Again, calculations show that this is not the case and we obtain a 
covariant version of the above integrability condition in solved form given 
by 

T _ AR) , AR) , AR) _ AR) _ AR) _ AR) _ AR) 

-'Coj; — ji212;33 ./1313;22 "T i2323;ll ./1323;12 ./1323;21 ./l213;23 ./l213;32 
+/l223;13 + /l223;31 • i^'^) 

This condition can be rewritten again in solved form as 

ICov = -Bi2123;3 + -Bi3i32;2 + -B23231;l i (58) 

( 

where Babcde '■= fab[cd-e] '^^ich amounts to the covariant derivatives of the 
Bianchi identities for the Rabcd involved in (]5^) . We can rewrite Icov [ P^ ] 
in a more concise form using bivectors, where we introduce the bivector- 
indices A := 12, S := 31, C := 23 . The integrability condition Icov can 
then be expressed as 

— fiR) I fiR) I f{R) \')(A^) I f(^) I A^) \ 

^Cov — ./^;(33) -r JbB;{22) -I CC_;ill) ^^■'AC;(13) -I BC;{12) ■'AB;{23)^ 
x,y=lX,Y=l 

Our new covariantly prolonged system is then given by 

n - f^^^ 

^ ~ J abed 

n - f 

~ Jabcd;e 

= Icov ■ (60) 

It is r{n^2^) = 6 + 19 = 25, r(7^^^^) = 6 + 19 + 38 = 63 as well as r{M^^^) 
= 19,r(A^^^^) = 38 so that m''^^ coincides with the previous symbol for 

(2) 

the prolongation involving partial derivatives. Again, it is dim(7e^ = 
dim(7^2^^) = 55. We find that the system ( pO| ) also consists of a system in 
involution and in addition to this also obeys general covariance. Therefore, 
we prefer ( |60[ ) as a prolongation to a second-order system in involution 
for the Riemann-Lanczos problem in 3 dimensions and we state the final 
proposition: 

Theorem 7 Proposition 

The Riemann-Lanczos problem in 3 dimensions is not in involution. Its 
reduced characters {s'q, s'^, S2, s'^) are (14,8,7,3). The differential gauge con- 
ditions do not change the non-involutivity and they modify the reduced char- 
acters {sq, s'l, S2, s'^) only slightly to (17,8,7,0). 

The 3-dimensional Riemann-Lanczos problem becomes involutive after 
just one prolongation which is obtained by adding either of the integrability 
conditions ( |5^ or (58) and now its Cartan characters {a'^\a^2\'^'2^) f^''" 
762^ are (18,9,2). 
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D. A Singular Solution for the 3-dimensional Reduced Godel 
Spacetime 

Lastly, we give a singular solution for the unprolonged 3-dimensional Godel 
spacetime. In 3 dimensions, the reduced Godel spacetime can be charac- 
terised by the following line element 

ds^ = a^{dt^ - dx'^ + ^e^^dy^ + 2e'dtdy) . (61) 

Here, we use the notation := t,x^ := x^x"^ := y for convenience. The 
Riemann-Lanczos equations together with the 3 differential gauge conditions 
are given explicitely in and are very similar to the equations in 4 dimen- 
sions. Here, we want to find a solution which only depends on x leading to 
the Ansatz 



Ltxy 


= Cie- , 


Ltyx 


= Cae- , 


Ltxt 


= C3, 


Lxyy 


= Cre^^ 


Ltxx 


= Qe- , 


Ltyt ■■ 


= Ce , 


Ltyy 


= Cse- , 


Lxyx 


= Cs , 



where Ci, ...Cg are arbitrary constants. When inserting this into the 6 inde- 
pendent Riemann Lanczos equations and the 3 differential gauge conditions, 
we obtain a solution which does not satisfy Lty^.^ = @ but when for in- 
stance choosing C5 = Ce = 0, we obtain the solution 

r, — -olpi^)'^ T, — oLpX 

^txy — g , J^tyx — g , 

T _ _g? T _ 3^ 2x 

^txt — g , ^xyy — 16 ^ 

Its singular solution manifold possesses 3-dimensional tangent spaces which 
are spanned by 



U 



d_ 

dt ' 

d a'^ ^. d a? ^ d 2)0? 2x ^ 



^ = 7^ + -tc'^'ttt h —e-"— h — — e 



dx 8 dLtxy 8 dLtyx 8 dLxyy 
1 ^,5a 1 9^ 3a , d 1 ^ ,5a 1 n,. 3a , d 



3 '16 2 16' dPtxyx 3 '16 2 16' dPtyxx 



3 Or d 

+-:-e 



4(Z dPxyyx 

We conclude that for this solution manifold sq = Sq = 3 while all other 
characters vanish. 



Conclusion 

In 2 dimensions, the Riemann-Lanczos problem is very simple and we showed 
that it is always in involution. The general solution was given for both 
possible choices of signatures, Lorentzian and Euclidean. 
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In 3 dimensions, a prolongation becomes necessary to make it a system 
in involution. An integrability condition based on the derivatives of the 
Bianchi identities occurs when we use Janet-Riquier theory and introduce 
the second-order partial derivatives Sahcde as new jet coordinates. A singular 
solution for the unprolonged problem for the reduced Godel spacetime is 
given. 
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Appendix A: Calculations for 2 Dimensions 

In this section we exhibit results and calculations which are necessary for 
the Riemann-Lanczos problem in 2 dimensions. First, i^a, • • • ,^6 are given 
by 



^3 = g^\2T\^ + r?2) + <7''(3rl2 + Tl^) + 5''r^2 , 
h = 9^\Ai + 3r?2) + ^''(rL + 2Tl^) + 5"rfi , 



(A.l) 



-55 



Pi2iig^\ + Pi22ig^\ + Pi2i2g^\ + Pi222g'^\ 



L^2^ 5- 



12103,1 



-^12115'"'^ 2 + -Pl22l5'^ 2 + -^12125^ 2 + -^12225'^ 2 — Ll21^: 
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12 



-22 



12103,2 



Ll22^A,l 
-^122^4,2 



The polar matrices when the differential gauge condition is included then 
look like 



equation 


dx^ 


dx'^ 


dLi2l 


dLi22 


dPl211 


dPl212 


dPl221 


dPl222 


"il212 
dLi2^.g 

K12I 
K\22 

2{U\dKi2i) 
2{U\dKi22) 
2{V\dKi2i) 
2{V\dKi22) 


"12121 
^5 

— -P12II 
— -P122I 

— C^1211 

— C^1221 

— V12II 

— ^1221 


"12122 

— -P1212 
— -P1222 
— f^l212 
— f^l222 

— V1212 

— ^1222 


2(ri2 — 1^22) 

-63 

1 












-mi+ri2) 

-64 


1 








9" 









-2 

g'' 







y2 



2 

g'' 

















where U forms a 1-dimensional integral element and span{U, V} form a 2- 
dimensional one. The linear multipliers Ai, ■ ■ ■ , Aq and Bi, ■ ■ ■ , Bq are given 
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by 



S7 



A^ = ^(yi - 5"^2), A, = -^g^^A^, 

1/1 1/2 

^1 - 2-^2(9 - -fjrg j, -02 - — ^ — , 

= i?2(53 - (rl^ - Tl^){g^^ - ^gii)), 54 = ^2(^4 + (rli + rf2)(5^2 _ ^^ii))^ 

= Be = ^{V^ - g^'B2), 

where 

^r = Ss = 2,- - - . 



Appendix B: Calculations for 3 Dimensions 

The quantity ctabcde for the Riemann-Lanczos problem in 3 dimensions is 
defined as 

O^abcde = ^ad,ei-^nbc + Lncb) + ^bc,ei^nad + L^da) ~ ^ac,e{^nbd + L^db) 
~^bd,ei^nac + Lnca) ■ 

We obtain 18 different otabcde for the 3-dimensional Riemann-Lanczos equa- 
tions when using the full polar matrix, which are given as follows Q 
(for e = 1,2,3): 

ai212e = -Rl212,e + "12126 + '^{^\2 + r22)-Pl21e — 2(rJ]^ -|- r^2)-fl22e 

+2r22-Pl31e + 2rf^P232e + 2rf2-pL23e — 4rf2-Pl32e 
ai313e = -Rl313,e + ai313e + '^^^■iPl21e + 2(r]^3 -|- r|3)Pi3ie — 2(rJ^ -|- rf3)Pi33e 

—2V\iP2ZZe — 4r^3Pi23e + 2r^3Pi32e 
a2323e = -R2323,e + a2323e — 2r33Pi22e — 2r22-Pl33e + 2(r23 -|- r33)P232e 

~2(r22 + r23)-P233e + 2r23i-'l23e + 2r23-Pl32e 
ai213e = -Rl213,e + ai213e + (1^13 + 2r23)-Pl21e — ri3-Pl22e + (1^12 + 2r23)Pl31e 

~ri2-Pl33e - ril-P232e — (Xll + 21^x2 " ri3)-Pl23e — (1^11 + 2rf3 — rf2)-Pl32e 
ai223e = -Rl223,e + ai223e — 1^23-^1216 + + r23)-fl226 + 1^22-^1316 — 1^22-^1336 

"(^12 + 2r'^3)P2326 + ri2-P2336 — (1^12 + 2r22 + r23)-f'l23e 

+ (2r23 + 2r22 ~ rj2)-Pl326 
ai3236 = -Rl323,6 + ai3236 — r33-Pl2l6 + r33-Pl226 + 1^23-^1316 — {'^'^12 + r23)-Pl336 

+ri3-f2326 — (ri3 + 2rf2)-P2336 — (ri2 + 2r23 — r|3)Pi23e 

+ (rl3 + 2ri3+r23)Pi32e. (B.2) 

The polar matrix of H(E'^)x for the 3-dimensional Riemann-Lanczos prob- 
lem can be written as 

/Moi Mo2Mo3\ 
\0 Ml ) ' 
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Mqi then is (for dx^ , dx"^ , dx^) the 6x3-matrix 



^ -Rl212,l + ai2121 
-Rl313,l + Q;13131 
-^2323,1 + a23231 
^1213,1 + ai2131 
-Rl223,l + ai2231 

\-Rl323,l + ai3231 



-Rl212,2 + «12122 
-Rl313,2 + Q!13132 
-^2323,2 + Q!23232 
-Rl213,2 + "12132 
-Rl223,2 + "12232 
-Rl323,2 + "13232 



-^1212,3 + "12123 \ 
-Rl313,3 + "13133 
-^2323,3 + "23233 
-Rl213,3 + "12133 
-^1223,3 + "12233 
-^1323,3 + "13233 / 



followed by the 8 columns (for the dLabc) forming the 6x8-matrix Mq2 which 
is 



2) -2{rli +r?2) 





-2rJ 
-r 



P 
13 „ 



2(r;3 



-2(rh +r3 
-2rl, 



-(2r}2 





■^(^23 ^ ^3; 
_r2 

-(rfs + 2rj 

^ 13 





-2r?i 
) -2(ri2 + ri3) 



2rf2) 



-(ri, 
-(-rL 



2r, 
■f 2r 



Z _|_ i-ij 



) 

J 

33) 



2rj3 - ri 
-ri, -r}. 



then followed by the 24 columns (for the dPabcd) producing the 6x24-matrix 



/0-2 200 00000 0000000 0\ 
00 -2 20000 0000000 
000000 00000-2 0200000 
-100 0-1 00000 0001001 
00-10 00010 0000200-1 

\0 000000 01000 1000010 0-2/ 



Ml consists of a 32x35-matrix which can be split into 



Ml := 



(M{P)l \ 
M{U) A^(;7) 
M{V) nIv) 

\mIz)qn\z) I 



M(P) is the 3x8- matrix 



M(P) := 



/ 


— -Pl211 


— A212 


— -P1213 \ 




— -fl221 


— -fl222 


— -fl223 




— -Pl311 


— -Pl312 


--Pl313 




— -P133I 


— -Pl332 


— ^1333 




— -f2321 


— -P2322 


— -f2323 




—^2331 


— -P2332 


—^2333 




-P123I 


— -P1232 


— -P1233 


V 


— -P132I 


— -F1322 


—^323 / 



M(U) is the 3x8- matrix 



M{U) := 



/ 


— ^^1211 


— ^^1212 


— ^^1213\ 




— t^l221 


— t^l222 


— t^l223 




— f^l311 


— f^l312 


— f^l313 




— t^l331 


— t^l332 


— f^l333 




— ^^2321 


— ^^2322 


— 1^2323 




— 1^2331 


— 1^2332 


— 1^2333 




— ^^1231 


— Ul232 


— ^^1233 


V 


— t^l321 


— t^l322 


— t^l323/ 



27 



N(U) then is the 8x24-matrix 

/[/I [/2J73 000000000000000000000\ 
000 C/i[/2[/3 000000000000000000 
000000 f/i[/2[/3 000000020000000 
000000000 f/iC/2c/3 000000000000 
000000000000 f/if72[/3 000000000 
000000000000000 t/iC/2C/3 000000 
000000000000000000 f/iC/2[/3 000 

Voooooooooooooooooooooc/ic/2 C/'V 

M(V) and M(Z) arc the same as M(U) only for the Vessiot vector fields V and Z. 
In the same way N(V) and N(Z) are the same as N(U). 

Next, we give the linear multipliers creating a linear combination for the reduced 
polar matrix. 







A2 = 




^3 = 




A4 = S2B1 , 




A,= 


yl i>l , 


^6 = 


y^<52l54 D 


Bi = Bi{arbitrary) , 




B2 = 


yrBi , 


B3 = 








B5 = 


V^l54 R 


Be = 


V^<54 R 


B7 = ^(|-2(54)i?i, 




Bs - 








Ci = -yrBi , 




C2 = 


R 




(7^<52 R 






C5 = 


R 


^6 = 


(7^54 R 




-2^U^Si)Bi,Cs^ 








where here 













and 






For the special choice Bi — U^V^V'^, we obtain a similar kind of identity as the 
one mentioned in (23) in |^ for the 4-dimensional case but not with exactly the 
same multipliers because we incorporated the cyclic conditions (|^) as well. 

Appendix C: A Coordinate Transformation for the 
3-dimensional Case 



The Bi in (|51|) are the parts which do not involve any term of the form Vabc3 and 
are given by 

Bl — a2ll4211 + 122^212 ~ 111^331 ~ O'ViV VlTl 
B2 ~ a31^1311 + ^32^2322 — aill4331 ~ ai2l^2332 
Bz ~ a3iV2321 + a32V2322 " 121^2331 — 022^2332 

-B4 — a'i\V\2\\ + a32Vl212 — 011^1231 ~ 022^1232 + a2lVi3ii + 022^1312 

— ailVl321 — ai2Vl232 
-B5 = 0-Z\yV22\ + 032^4222 — 2a2lPl231 — 2022^^1232 + a2lVl321 + 022^1322 
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— 111^2321 — 112^2322 

Be = 2a3iVi32l + 2a32Vi322 — a2lVi33i — a22Vi332 — 031^1231 " 032^1232 

— 011^2331 ~ 012^2332 ■ 

The function /(K6ci,K6c2) for ^2332 = /(Vahci, K6c2) is given by 

2 2 

X. _ ^23 f^r / 031033 021033 N , ^033033 022033 X 

1^2332 — [l'l21l(, ) + >'l212l, ) 

O23O12 — O22O13 ai3 013023 O13 013023 

, y , .011033^ 031033 .012033^ O32O33 

+ V'l22l(, j + 1'1222(, j + V1311 

O13O23 O23 O13O23 023 

/O21O33 021033 ,022033 O23O33 ,011033 

( ) + Vi312( ) + T433l( 021) 

Ol3 O13 O13 O23 Oi3 

, Tr / \ I T/ /O11O33 O13O31 ,012033 O13O32, 

+ l'1332(,Ol2 — 022j + l'232l(, j + ''2322 j 

O23 O23 O23 O23 

, ,013021 , ,„02l033 O11O33 

+ v'233l(, Oiij + Vl23ll^ 031; + V1232 

O23 O23 023 

,r,022033 O12O33 , , O11O33 O33O21, 

(2 ) + Vl32l(2o3i ) 

O23 O13 — O32 O13 O23 

IT/ O12O33 O22O33 

+ Vl322(2032 ) . (C.2) 

Ol3 O23 



^ All monomials which are not in the set of principal monomials or in its extension 

by multiplicative variables correspond to parametric derivatives. 

^ Note that if a system of partial differential equations is translated into an exterior 

differential system in the appropriate way, then we have the correspondence 

Sk = OL'q^ as long as it is in involution. We refer the reader to Q for further 

details concerning involution. 
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